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Abstract 

In  this  note,,  we  update  ^earlier  results  on  how  various  types  of  Coxian  distributions 
relate  to  each  other.  Some  more  recent  uses  of  such  distributions  are  discussed,  with 
an  emphasis  on  their  connection  to  generalized  hyperexponential  distributions.  We  pay 
particular  attention,  to  the  class  of  densities  having  Laplace  transforms  with  only  real 
negative  zeros  and  poles  and  their  relationship  with  the  GH  class. 
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1 .  Background 

A  number  of  recent  papers  (e  g.,  Ott,  1987,  and  Shanthikumar,  1985)  have  used 
distribution  functions  related  to  the  classical  exponential  staging  formulation  introduced 
by  Erlang  (see  Brockmeyer  et  al.,  1948),  with  primary  modifications  over  the  years  by 
Jensen  (1954)  and  Cox  (1955).  Much  of  the  current  popularity  of  such  distributions  is 
due  to  the  work  of  Neuts  and  colleagues  on  the  so-called  phase-type  distributions  (see, 
e.g.,  Neuts,  1975a, b  and  1981),  exploiting  relationships  to  the  theory  of  Markov  chains 
and  putting  the  theory  to  effective  computational  use  in  a  large  variety  of  stochastic 
models.  Some  recent  work  by  these  authors,  together  with  a  number  of  different 
coauthors,  has  focused  on  the  generalized  mixed  exponential  form  (called  GH)  of 
Coxian  distribution  (see  Botta  &  Harris,  1986,  Botta,  Harris  &  Marchal,  1987,  Harris  & 
Sykes,  1987).  These  are  linear,  but  not  necessarily  convex  combinations  of  negative 
exponential  densities.  Botta  and  Harris  (1986)  showed,  critically,  that  the  GH  class  is 
dense  in  the  set  of  all  CDFs  relative  to  an  appropriate  metric.  Denseness  is  likewise  a 
property  of  both  the  PH  and  Coxian  classes  (the  latter  henceforth  called  Rn). 

For  purposes  of  clarification,  we  note  the  following  relationships  between  the  major 
classes  of  cumulative  distribution  functions  (CDFs).  For  the  most  part,  all  of  this 
stems  from  Erlang's  early  idea  of  modeling  a  duration  or  lifetime  as  a  sum  of  (k)  in¬ 
dependent  and  identical  exponential  stages.  Much  later,  Jensen  (1954)  generalized 
Erlang’s  device  to  allow  the  stages  to  have  non-identical  CDFs,  and  indeed  recognized 
the  natural  connection  between  Erlang’s  method  of  stages  and  absorption-time  distribu¬ 
tions  for  finite  Markov  chains.  (Note  that  when  there  is  possible  exit  from  any  stage, 
an  Erlang  becomes  a  random  sum  of  exponentials.)  An  important  milestone  in  this  de¬ 
velopment  came  very  shortly  thereafter  in  work  by  Cox  (1955),  who  generalized  to 
cover  all  distributions  with  rational  Laplace  transforms,  by  involving  a  more  compli¬ 
cated  stage-to-stage  movement  possibly  using  complex  transition  probabilities  and  even 
an  infinite  number  of  steps.  More  precisely,  the  sequential  flow  corresponding  to  a 
Coxian  distribution  can  have  negative  branching  probabilities  and  complex  scale  pa¬ 
rameters  (with  negative  real  parts).  While  such  stages  may  not  have  a  physical  reality, 
the  resulting  CDF  can  well  be  legitimate.  The  class  Rn  contains  the  class  of  all 
phase-type  distributions,  which  in  turn  contains  all  generalized  Erlangs. 

Recently,  Shanthikumar  (1985)  worked  with  two  classes  of  functions,  which  he 
called  generalized  and  bilateral  phase  types.  The  CDF  of  the  generalized  phase  type 
(GPH)  is  created  from  an  infinite  mixing  on  the  number  of  convolutions.  The  bilateral 


phase  type  is  defined  over  the  entire  real  line  in  an  analogous  fashion  using  Erlang 
mixing.  It  follows  that  a  GPH  distribution  is  an  ordinary  PH  whenever  the  mixing  dis¬ 
tribution  has  finite  support  or  has  an  infinite  phase-type  representation. 

Ott  (1987)  has  renamed  the  GPH  distributions  as  almost  phase  types  (APH),  and 
has  developed  a  relationship  between  the  Wiener-Hopf  factorization  for  the  G/G/l 
queue  and  infinite  matrix  (possibly  phase-type)  representations  of  at  least  one  of  the 
interarrival  and  service  CDFs.  An  iterative  numerical  procedure  is  formulated  for  solv¬ 
ing  the  special  cases  of  the  GPH/G/1  and  G/GPH/1  queues.  In  these  cases,  the  algo¬ 
rithm  works  from  a  finite  mixture  of  Erlangs,  even  though  this  might  be  just  one  ele¬ 
ment  in  an  infinite  sequence.  Since  the  GPH  class  is  dense  in  the  set  of  all  CDFs, 
one  can  theoretically  always  find  such  a  close  approximant. 

In  another  recent  paper,  Swensen  (1986)  dealt  with  a  variation  on  the  GI/M/c 
model  in  which  customers  renege  when  their  waiting  times  exceed  a  fixed  value.  To 
make  the  problem  as  general  as  possible,  the  CDF  for  interarrival  times  was  assumed 
to  belong  to  a  general  class,  said  to  be  Coxian.  However,  by  restricting  the  mixing 
probabilities  and  means  of  the  exponential  stages  to  positive  values,  Swensen’s  class  of 
CDFs  is,  in  fact,  identical  to  the  class  of  finite  mixtures  of  generalized  Erlang  distribu¬ 
tions.  While  this  restricted  class  is  a  subset  of  the  PH  class  of  distributions,  the  larger 
class  of  Coxian  distributions  actually  contains  the  PH  class,  contrary  to  Swensen’s  state¬ 
ment  that  the  Coxian  distributions  are  a  subset  of  the  phase-type  CDFs. 

Although  these  observations  do  not  affect  any  of  the  results  in  the  paper,  it  is  im¬ 
portant  to  note  that  when  the  author  speaks  of  Coxian  distributions,  he  is  actually  con¬ 
sidering  only  a  very  restricted  subset  of  those  distributions.  Confusion  may  result  if 
this  limited  definition  of  Coxian  distributions  is  not  borne  in  mind  when  interpreting 
Swensen’s  results. 

The  newest  class  added  to  the  list  has  been  offered  by  Sumita  and  Masuda  (1987). 
This  is  the  class  (called  D+)  of  probability  distribution  functions  which  have  Laplace- 
Stieltjes  transforms  with  only  real  negative  zeros  and  poles.  These  densities  are  the 
newest  and  therefore  the  least  well  known,  and  because  their  possible  use  raises  some 
interesting  questions,  we  devote  the  next  section  of  the  note  to  a  discussion  of  the  po¬ 
sition  of  Q+  relative  to  the  other  types  of  phase  classes. 


2.  The  Density  Class  Cl+ 

The  formal  definition  of  f2+  is  the  class  of  densities  admitting  Laplace  transforms 


n  (1  +  s/ti  ) 
i=i  1 

1  .  "  “  1  » 

n  (us/Bi) 

j=  i  1 


0  <  m  <  n  <  oo. 


^i'®j  >  0  for  al1  ij- 

(We  assume  that  the  numerator  is  1  for  m  =  0.] 

Clearly,  all  such  transforms  are  rational  and  thus  Cl  c  Rn  .  Since  all  {0  j }  are  positive, 
it  appears  at  first  glance  that  f*(s)  should  also  be  generalized  hyperexponential  when¬ 
ever  the  {0  j }  are  district.  However,  this  is  not  the  case. 

Sumita  and  Masuda  (1987)  have  actually  given  the  skeleton  of  a  counterexample 
pdf  which  is  in  GH  but  not  to+.  Consider  the  transform 

f  *  ( si  —  ^  2s2+10s+l  3 

13  (s+l)(s+2)(s+3) 

This  is  easily  inverted  to  be  the  GH  density 


f(t)  =_j!(e")--l(2e-’‘)  +  -l(3e-n 


However,  we  see  that  f(t)  is  not  in  since  the  roots  of  the  numerator  of  f*(s)  are 
complex.  This  particular  example  arose  as  a  mixture  of  two  densities,  thus  showing 
that  G+  is  not  closed  under  mixing. 

Clearly,  this  is  not  a  pathological  counterexample,  for  there  are  many  examples  of 
GH  distributions  which  have  similar  transform  properties.  Another  illustration  is  Exam¬ 
ple  2.2.2  (page  124)  of  Botta,  Harris  and  Marchal  (1987),  namely, 

f(t)  =  4(e-‘)-6(2cr2l)+3(3<r3t). 

The  numerator  of  its  transform  is  the  polymomiai 

fn(s)  =  s  -s+6  , 

which  also  has  complex  roots.  Furthermore,  similar  complex  arithmetic  issues  create 
examples  of  distributions  in  the  Neuts  phase-type  class  PH  which  are  also  not  in  Q+. 
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Some  of  the  major  results  of  Sumita  and  Masuda  (1987)  have  interesting  applica¬ 
tion  to  the  theory  of  Coxian  and  generalized  hyperexponential  distributions.  A  prime 
example  of  these  is  their  Theorem  1.2,  which  provides  a  simple  sufficient  condition  for 
f  e  GH  when  f  is  also  in  fl+.  We  have  been  able  to  extend  their  argument  to  derive  a 
sufficient  condition  for  f  e  GH  independent  of  whether  it  is  in  . 

3.  Sufficient  Conditions 

The  condition  offered  by  Sumita  and  Masuda  requires  that  there  exists  an  indexing 
of  the  {0j}  and  the  {T)j}  (assuming  tj ,  <  r|2  <  .  .  .  <  and  0,  <  02  <  .  .  .  <  0n)  in 

which  0 ;  <  T]  j  for  1  <  i  <  m.  An  informal  proof  of  this  result  follows  from  writing 
the  transform  as  the  product  of  quotients 

1  ±  s/n,  •  •  •  l_±_sZllm  _ 1 _  •  •  •  _ 1 _ • 

1  +  s/0,  1  +  s/02  1  +  s/0m  1  +  s/0  m+l  1  +  s/0n 

The  inverse  Stieltjes  transform  of  each  of  the  m  factors  is  of  the  form 

F(t)  =  +  (1  -  £j)(l  -  e-0il  )  . 

T)  i  "H  ' 


The  requirement  that  0  j  <  t]  4  yields  a  mixture  of  an  atom  at  the  origin  and  an  expo¬ 
nential. 

Clearly,  each  of  the  last  n-m  terms  corresponds  to  an  exponential.  Since  each  of 
the  n  factors  thus  corresponds  to  a  legitimate  probability  distribution,  the  convolution 
of  the  corresponding  functions  will  yield  a  true  distribution  (without  an  atom  since 
there  is  at  least  one  purely  exponential  term  in  the  convolution). 

This  suggests  the  following  extension  to  the  case  where  the  numerator  polynomial 
of  the  transform  can  have  complex  roots  occurring  in  conjugate  pairs. 

Theorem  Suppose  a  rational  transform  has  the  form 

m 

n  (i+s/0,  ) 

f * (s)  =  — -  (0  <  m  <  n), 

n 

n  (l+s/e. ) 

j=i  1 
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i.'iVi 
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where  the  {0j}  are  real,  positive  and  arranged  in  ascending  order,  and  the  {r|  i)  are 
either  real  and  positive  or  occur  in  complex  conjugate  pairs  with  positive  real  parts. 
Suppose  (without  loss  of  generality)  that 

Re  Cn,)  <  Re  (“n2)  <  .  .  .  <  Re  (r)m)  [where  equality  holds  only  in  the  case  of 

complex  conjugates] 

and  that  for  i  =  1,  2,  .  .  .  ,  m,  0 ;  <  q  ;  when  -q  s  is  real  and  (9  i  +  6  j+i )/2  <  Re(n  j) 
when  (rij,  T|  i+1)  are  a  complex  conjugate  pair.  Then  the  inverse  transform  of  f*(s)  is  a 
probability  distribution. 

Proof:  The  idea  of  the  proof  is  the  same  as  before,  except  that  we  now  can  have  fac¬ 
tors  of  the  form 

)  (l+-\) 

'  a+zb  '  '  a-/b ' 

(1+  s/rj  (i+  s/r2) 

Each  of  such  factors  can  be  expressed  as 


2a-r.  -r2  r.r. 

- 1 - 2 -  S  +  1  _  1  2 

2  u2  o2  K2 

a  +b _ a  +b 

(1  +  s/r,  )(1  +  s/r2) 


2  i2 

a  +b 


Expanding  the  second  term  in  partial  fractions  yields 


r,2-2ar,+a2+b2 


Since  r2  <  r, ,  e  >  e  and  thus  the  term  in  the  braces  will  be  nondecreasing  if 


(r2-a)2  <  (r,-a)2 


Now,  if  a  <  r2  >  r,,  this  inequality  is  satisfied; 
if  r2  >  r,  >  a,  the  inequality  is  violated;  and 
if  r2  >  a  >  rlt  the  inequality  is  satisfied  as  long 

as  a  >  (r,  +  r2)/2.  Therefore  a  >  (r,  +  r2)/2  assures  that  F(t)  is  nondecreasing  and,  since 
F(oo)  =  1,  it  is  a  legitimate  CDF. 

Real  terms  are  tested  as  before,  so  that  every  factor  corresponds  to  a  probability 
function,  and  the  resulting  convolution  is  also  probability  function  (again  with  no  atom 
at  the  origin).  Consider  the  following  example.  First,  let 


f*(s)  = 


(>■ 


2+2/ 


'  v  2+2  /  ' 


(l+s/l)(l+s/2)(l+s/3) 


The  inverse  of 


(u-J— )(u-I_) 

v  2 +/2/v 


2-/2 


(  1+s/l)  (  l+s/2) 


is 


F,(t)  -U5  (l-e"‘  )  -I  0-e-2'  ). 

4  4  2 

This  is  indeed  a  CDF  since 

61+62  =  i.  <  Re(rj,)  =  2. 

2  2 

We  also  see  that  the  inverse  of  1/(1  +s/3)  is  F2(t)  =  1  -  e-3t  ,  and  it  follows  that  F(t)  is 
a  CDF. 


-  7  - 


As  a  second  illustration,  let  F2(t)  = 

1  -  e-t  +  e-21 

-  e->‘  , 

with 

f(s)  =  J_ 

-  J_  +  3 

1+s 

2+s  3+s 

=  6 

+  6s  +  2s2 

0  + 

s)(2+s)(3+s) 

f, 

s 

¥  i  *  _s 

A  (3— s/3  i)/2 

.  V_ 

(3+\/3  i)/2 

(]+s)(2+s)(3+s) 

Taking 

3+^3  i 

1.  r2  =  2, 

’  2 

■  i  ^2  -  "n  i » fi  - 

we  have 

LIU  3  = 

Re  (tj  , ) 

2  2 

so  that  the  condition  is  satisfied  and  F(t)  is  therefore  a  CDF. 
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